Results of the search for inspiraling compact star binaries 
from TAMA300's observation in 2000-2004 



Tomomi Akutsu, 1 Tomotada Akutsu, 2 Masaki Ando, 3 Koji Arai, 4 Akito Araya, 5 Hideki Asada, 6 Youichi Aso, 3 Mark 
A. Barton, 1 Peter Beyersdorf, 4 Youhei Fujiki, 7 Masa-Katsu Fujimoto, 4 Ryuichi Fujita, 8 Mitsuhiro Fukushima, 4 

Toshifumi Futamase, 9 Yusaku Hamuro, 7 Tomiyoshi Haruyama, 10 Hideaki Hayakawa, 1 Kazuhiro Hayama, 2: LJ 
Gerhard Heinzel, 11 Gen'ichi Horikoshi, 10 'Q Hideo Iguchi, 12 Yukiyoshi Iida, 3 Kunihito Ioka, 13 Hideki Ishitsuka, 1 
Norihiko Kamikubota, 10 Nobuyuki Kanda, 14 Takaharu Kaneyama, 7 Yoshikazu Karasawa, 9 Kunihiko Kasahara, 1 
Taketoshi Kasai, 6 Mayu Katsuki, 14 Keita Kawabe, 15 Mari Kawamura, 16 Seiji Kawamura, 4 Nobuki Kawashima, 17 
Fumiko Kawazoe, 18 Yasufumi Kojima, 19 Keiko Kokeyama, 18 Kazuhiro Kondo, 1 Yoshihide Kozai, 4 Hideaki Kudoh, 20 
Kazuaki Kuroda, 1 Takashi Kuwabara, 7 Namio Matsuda, 21 Norikatsu Mio, 22 Kazuyuki Miura, 23 Osamu Miyakawa, 24 
Shokcn Miyama, 4 Shinji Miyoki, 1 Hiromi Mizusawa, 7 Shigenori Moriwaki, 22 Mitsuru Musha, 25 Shigeo Nagano, 26 

Yoshitaka Nagayama, 14 Ken'ichi Nakagawa, 25 Takashi Nakamura, 16 Hiroyuki Nakano, 14 'Q Ken-ichi Nakao, 14 
Yuhiko Nishi, 3 Kenji Numata, 27 Yujiro Ogawa, 10 Masatake Ohashi, 1 Naoko Ohishi, 4 Akira Okutomi, 1 Ken-ichi 
Oohara, 7 Shigemi Otsuka, 3 Norichika Sago, 8 '0 Yoshio Saito, 10 Shihori Sakata, 18 Misao Sasaki, 28 Kouichi 
Sato, 29 Nobuaki Sato, 10 Shuichi Sato, 4 Youhei Sato, 25 Hidetsugu Seki, 3 Aya Sekido, 30 Naoki Seto, 31 Masaru 

Shibata, 32 Hisaaki Shinkai, 33 Takakazu Shintomi, 10 Kenji Soida, 3 Kentaro Somiya, Toshikazu Suzuki, 10 
Hideyuki Tagoshi, 8 Hirotaka Takahashi, 34 Ryutaro Takahashi, 4 Akiteru Takamori, 5 Shuzo Takemoto, 16 Kohei 
Takeno, 22 Takahiro Tanaka, 16 Keisuke Taniguchi, 35 Shinsuke Taniguchi, 3 Toru Tanji, 22 Daisuke Tatsumi, 4 C.T. 

Taylor, 1 Souichi Telada, 36 Kuniharu Tochikubo, 3 Masao Tokunari, 1 Takayuki Tomaru, 10 Kimio Tsubono, 3 
Nobuhiro Tsuda, 29 Yoshiki Tsunesada, 4 Takashi Uchiyama, 1 Akitoshi Ueda, 4 Ken-ichi Ueda, 25 Fumihiko Usui, 37 
Koichi Waseda, 4 Yuko Watanabe, 23 Hiromi Yakura, 23 Akira Yamamoto, 10 Kazuhiro Yamamoto, 1 Toshitaka 
Yamazaki, 4 Yuriko Yanagi, 18 Tatsuo Yoda, 3 Jun'ichi Yokoyama, 38 Tatsuru Yoshida, 9 and Zong-Hong Zhu 4 

(The TAMA Collaboration) 

institute for Cosmic Ray Research, The University of Tokyo, Kashiwa, Chiba 277-8582, Japan 
Department of Astronomy, The University of Tokyo, Bunkyo-ku, Tokyo 113-0033, Japan 
Department of Physics, The University of Tokyo, Bunkyo-ku, Tokyo 113-0033, Japan 
^National Astronomical Observatory, Mitaka, Tokyo 181-8588, Japan 
5 Earthquake Research Institute, The University of Tokyo, Bunkyo-ku, Tokyo 113-0032, Japan 
6 Faculty of Science and Technology, Hirosaki University, Hirosaki, Aomori 036-8561, Japan 
7 Faculty of Science, Niigata University, Nngata, Niigata 950-2102, Japan 
8 Graduate School of Science, Osaka University, Toyonaka, Osaka 560-0043, Japan 
9 Graduate School of Science, Tohoku University, Sendai, Miyagi 980-8578, Japan 
10 High Energy Accelerator Research Organization, Tsukuba, Ibaraki 305-0801, Japan 
11 Max- Planck- Institut fur Gravitationsphysik, Callinstrasse 38, D-30167 Hannover, Germany 

12 Tokyo Institute of Technology, Meguro-ku, Tokyo 152-8551, Japan 
Physics Department, Pennsylvania State University, University Park, Pennsylvania 16802, USA 
•* Graduate School of Science, Osaka City University, Sumiyoshi-ku, Osaka 558-8585, Japan 
15 LIGO Hanford Observatory, Richland, Washington 99352, USA 
lb Faculty of Science, Kyoto University, Sakyo-ku, Kyoto 606-8502, Japan 
11 Kmki University, Higashi- Osaka, Osaka 577-8502, Japan 
1 Ochanomizu University, Bunkyo-ku, Tokyo 112-8610, Japan 
10 Department of Physics, Hiroshima University, Higashi- Hiroshima, Hiroshima 739-8526, Japan 
20 Theoretical Astrophysics Group, Department of Physics, 
The University of Tokyo, Bunkyo-ku, Tokyo 113-0033, Japan 
21 Tokyo Denki University, Chiyoda-ku, Tokyo 101-8457, Japan 
22 Department of Advanced Materials Science, The University of Tokyo, Kashiwa, Chiba 277-8561, Japan 
23 Department of Physics, Miyagi University of Education, Aoba Aramaki, Sendai 980-0845, Japan 
'^California Institute of Technology, Pasadena, California 91125, USA 
25 Institute for Laser Science, University of Electro -Communications, Chofugaoka, Chofu, Tokyo 182-8585, Japan 
20 National Institute of Information and Communications Technology, Koganei, Tokyo 184-8795, Japan 
27 NASA Goddard Space Flight Center, Greenbelt, Maryland 20771, USA 
28 Yukawa Institute for Theoretical Physics, Kyoto University, Sakyo-ku, Kyoto 606-8502, Japan 
29 Precision Engineering Division, Faculty of Engineering, 
Tokai University, Hiratsuka, Kanagawa 259-1292, Japan 
30 Waseda University, Shinjyuku-ku, Tokyo 169-8555, Japan 
Theoretical Astrophysics, California Institute of Technology, Pasadena, California 91125, USA 
32 Graduate School of Arts and Sciences, The University of Tokyo, Meguro-ku, Tokyo 153-8902, Japan 
33 Department of Information Science, Osaka Institute of Technology, Hirakata, Osaka 573-0196, Japan 



2 



* Max-Planck-Institut fur Gravitationsphysik, Albert- Einstein-Institut, 
Am Muhlenberg 1, D-14476 Golm bei Potsdam, Germany 
Department of Physics, University of Illinois at Urbana- Champaign, Urbana, Illinois 61801-3080, USA 
16 National Institute of Advanced Industrial Science and Technology, Tsukuba, Ibaraki 305-8563, Japan 
37 ISAS/JAXA, Sagamihara, Kanagawa 229-8510, Japan 
38 Research Center for the Early Universe(RESCEU), 
Graduate School of Science, The University of Tokyo, Tokyo 113-0033, Japan 

We analyze the data of TAMA300 detector to search for gravitational waves from inspiraling 
compact star binaries with masses of the component stars in the range 1 — 3Mq. In this analysis, 
2705 hours of data, taken during the years 2000-2004, are used for the event search. We combine the 
results of different observation runs, and obtained a single upper limit on the rate of the coalescence 
of compact binaries in our Galaxy of 20 per year at a 90% confidence level. In this upper limit, the 
effect of various systematic errors such like the uncertainty of the background estimation and the 
calibration of the detector's sensitivity are included. 

PACS numbers: 95.85.Sz,04.80.Nn,07.05.Kf,95.55.Ym 



I. INTRODUCTION 

Several laser interferometric gravitational wave detec- 
tors of the first generation have been operated to de- 
tect gravitational wave. These include GEO[l|, LIGOQ, 
TAMA3OO0, and VIRGO 0. These detectors have been 
improved their sensitivity very rapidly in the past sev- 
eral years. The direct detection of gravitational waves is 
important not only because it will become a new astro- 
nomical tool to observe our universe, but also because 
it will become a new tool to test general relativity and 
other gravity theories in a strong gravity field region. 

In this paper, we present results of the data analysis of 
the TAMA300 detector to search for gravitational waves 
produced by inspiraling compact star binaries, comprised 
of non-spinning neutron stars and/or black holes. Inspi- 
raling compact binaries are considered to be one of the 
most promising sources for ground based laser interfer- 
ometers. TAMA300 has been performed nine observation 
runs of the detector since 1999. The total amount of data 
is more than 3000 hours. Given such large amount of 
data, it is very interesting to analyze the data to search 
for candidate gravitational wave events and to set an up- 
per limit to the event rate. 

In the past, there were several works which searched 
for inspiraling compact binaries using laser interferome- 
ter data. Allen et al. @ analyzed LIGO-40m's data in 
the mass range 1 — 3Mq, and obtained an upper limit of 
0.5 [1/hour] on the Galactic event rate. The data from 
"Data Taking 2" (DT2) of TAMA300 in 1999 was ana- 
lyzed in the mass range 0.3 — 10M© |6j, and an upper 
limit of 0.59 [1/hour] on the event rate with signal-to- 
noise ratio greater than 7.2 was obtained. TAMA300's 
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DT6 data and LISM-20m's data taken in 2001 were an- 
alyzed to search for coincident signals, and obtained an 
upper limit of 0.046 [1/hour] on the nearby event rate 
within lkpc from the Earth @|. Abbott et al. Q an- 
alyzed LIGO's "1st Scientific run" (SI) data taken in 
2002, and obtained an upper limit of 1.7xl0 2 per year 
per Milky Way Equivalent Galaxy (MWEG) in the mass 
range 1 — 3M©. Abbott et al. analyzed LIGO S2 data 
taken in 2003, and obtained an upper limit of 47 per year 
per MWEG in the mass range 1 — 3M© |9(, 63 per yea r 
per Milky Way halo in the mass range 0.2 — 1M© |1CJ| . 
LIGO's S2 data was also analyzed to search for binary 
black hole inspirals in the mass range 3 — 20M©, and 
an upper limit of 37 per year per MWEG was obtained 
H3. LIGO's S2 data and TAMA300's DT8 data were 
analyzed to search for coincident signals and an upper 
limit of 49 per year per MWEG was obtained [l^ . In 
all of above cases, there were no signals that could be 
identified as gravitational waves. 

In this paper, we analyze the data from DT4, DT5, 
DT6, DT8, and DT9 of TAMA300. A part of DT6 data 
which was coincident with LISM was already analyzed in 
[7|. The initial results of the analysis of DT8 data was 
reported in Rcf. |H3. A part of DT8 data which was 
coincident with LIGO S2 was already analyzed in [l2l ]. 
In this paper, we analyze these data again together with 
the other data in a unified way. Until DT6 observation, 
TAMA300 was the only large scale laser interferometer 
which was operated. Thus, it is important to analyze 
such data to search for possible gravitational wave sig- 
nals. Further, in order to take advantage of the long 
length of data from DT6, DT8 and DT9, we combine the 
results from DT6, DT8 and DT9 data and obtain a sin- 
gle upper limit on the rate of the coalescence of compact 
binaries in our Galaxy. We also evaluate the systematic 
errors caused by the uncertainty of the calibration and 
the background trigger rate. Other errors such like the 
uncertainty of the distribution model of sources, and the 
uncertainty of the theoretical templates are also evalu- 
ated. These systematic errors are taken into account to 
evaluate the upper limit. 
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This paper is organized as follows. In section II, we 
overview the detector and the data we analyze. In section 
III, the analysis method is presented. In section IV, the 
results of the analysis are presented. In section V, the 
evaluation of the detection probability of the Galactic 
signals and the upper limit to the event rate are shown. 
In section VI, we evaluate the errors due to various error 
sources, and its effect to the upper limit. In section VII, 
we summarize the results and present the conclusion. 



II. DATA FROM THE TAMA300 DETECTOR 

TAMA300 is a Fabry-Perot-Michelson interferometer 
with baseline length of 300m located in Mitaka, Tokyo 
(35°40'N, 139°32'E). The history of the observation run 
of TAMA300 is listed in Tabled Until DT6, the detector 
was operated without the power recycling system. Af- 
ter DT6, the power recycling system was installed. The 
main signal of detector is recorded with a 20 kHz, 16 
bit data-acquisition system. There are more than 150 
signals which monitor the condition of the detector, and 
the environment of detectors. During the operation, the 
mirrors of the detector are shaken by a 625 Hz sinusoidal 
signal in order to calibrate the detector sensitivity con- 
tinuously. 

We use DT4, DT5, DT6, DT8, and DT9 data of 
TAMA300. The observations of TAMA300 were inter- 
rupted by the unlock of the detector. They were some- 
times suspended manually for maintenance. By removing 
such dead time, the total length of data available for the 
data analysis is 3032 hours. Among them, we do not 
use the first 6.5 minutes of data just after the detector 
recovers from the dead time, because such data often con- 
tain signals due to the excitation of the violin modes of 
pendulum wires, and/or other signals caused by distur- 
bance during the dead time. The data from the detector 
are transfered into the strain equivalent data by apply- 
ing the transfer function. The fluctuation of the transfer 
function at each time is determined by computing the 
optical gain. We do not use the data if the value of the 
optical gain deviates from the average value significantly. 
The total amount of data remained after removing such 
bad quality parts are 2705 hours. We analyze these data 
to search for gravitational wave events. However, we do 
not use DT4 and DT5 data to set the upper limit for the 
event rate, because the length of data from these runs are 
much shorter than DT6-8-9, and because quality of data 
of these runs are much worse than those of DT6-8-9. The 
total amount of data used for setting the upper limit is 
2462.8 hours. 



III. ANALYSIS METHOD 

The standard method to search for gravitational wave 
signals with known wave forms in noisy data is the 
matched filtering method, in which we search for the 



best matched parameters of the theoretical wave form 
by cross-correlating the data with the theoretical wave 
form. As the theoretical wave form, we use the non- 
spinning, restricted post-Newtonian (PN) wave form in 
which the phase is given to high post-Newtonian order, 
but only the leading quadrupole term is contained in the 
amplitude. We use the phase formula derived from 2.5 
PN approximation. Although the current best formula 
by the PN approximation is the 3.5 PN formula [Tj|, the 
error due to the use of 2.5 PN formula instead of 3.5PN 
formula is small for binaries of mass considered in this 
paper (see Section IV). On the other hand, the PN ap- 
proximation itself may contain errors due mainly to the 
relativistic effects in the region when the orbital radius is 
the same order as the gravitational radius of stars. These 
effects will be incorporated in the systematic error to the 
detection probability of signals. 

The basic formula of the matched filtering method is 
given by 

p = ^(s,h ) 2 + (s,K /2 ) 2 , (3.1) 

where 

M) = 2 Lw* (3 ' 2) 

and where a(/) and b(f) are the Fourier transformation 
of time sequential data, a(t) and b(t). The Fourier trans- 
formation is defined by 

/oo 
a{t)e 2 * lft dt. (3.3) 
-oc 

The function s(t) is the time sequential data from the 
detector. Two functions, ho and h^i^ are the tem- 
plates in the frequency domain, which are normalized 
as (ho, ho) = (hir/zjhir/s) = 1. The Fourier transforma- 
tion of them, ho(f) and h^^if), are computed by the 
stationary phase approximation. We thus have the or- 
thogonality of them, i.e., (ho,h w /2) = 0. The function 
S n (f) is the one-sided power spectrum density of noise 
of the detector. 

The parameters which describe the wave form are the 
time of coalescence, t c , the phase of wave at the coales- 
cence, 4> c , the total mass M = mi + mi and the non- 
dimensional reduce mass r\ = mi 777,2/ M 2 of the binary. 
We search for the parameters which give the maximum of 
p. In the formula (|3.1I) , the maximization over the phase 
is already taken analytically. The value of parameters, 
t c , M and 77, which maximize p are searched numerically. 

The data are divided into subsets of data with length 
52.4 seconds. Each subset of data has overlapping data 
with adjacent data for 4.0 seconds. Each subset of data 
is Fourier transformed, and the components more than 
5kHz are removed. The data are converted to the strain 
equivalent data s(/) by the transfer function. The power 
spectrum density of noise S n (f) is evaluated at neighbor 



4 



Period Typical strain noise Observed data Analyzed data 







[1/VHz] 


[hours] 


[hours] 


DTI 


6-7 Aug. 1999 


3 x 10 J 


11 




DT2 


17-20 Sept. 1999 


n ., ,-. — 90 

3 x 10 


31 




DT3 


20-23 April 2000 


H ., _ on 

1 x 10 


13 




DT4 


21 Aug.-4 Sept. 2000 


1 x 10 


154.9 


147.1 


DT5 


2-10 Mar. 2001 


1.7 x 10~ 20 


107.8 


95.26 


DT6 


1 Aug.-20 Sept. 2001 


5 x 10~ 21 


1049 


876.6 


DT7 


31 Aug.-2 Sept. 2002 




25 




DT8 


14 Feb.-14 April 2003 


3 x 10 -21 


1163 


1100 


DT9 


28 Nov. 2003 - 10 Jan. 2004 


1.5 x 10" 21 


556.9 


486.1 


Total 






3111.6 


2705 



(DT6,8,9 for upper limit) (2462.8) 



TABLE I: Observation history of TAMA300. For each observation run, the period of the observation, typical strain equivalent 
noise level around the most sensitive frequency region, the length of data observed, and the the length of data analyzed in this 
paper are shown. 



of each subset. Details of the method to evaluate S n (f) 
was described in §111. B. of Ref. Q. With the subset of 
data, we compute p. For each small time interval with 
length At c = 25.6msec, we search for t c , M and r\ which 
give the maximum of p. The value of p at all of t c can 
be computed automatically from the inverse FFT of the 
inner product, Eq. (|3.2fl . with respect to t c . The search 
for the best matched M and r\ is done by introducing the 
grid points in the two dimensional mass parameter space. 
The range of masses of each member star of binaries is set 
to 1 to 3 Mq. The grid separation length is determined 
so that the minimal match is less than 3% [l5|. The 
actual mass parameters we use for setting up the mass 
parameter space are those discussed in [l6|. We define 
a trigger by the local maximum of p in each small time 
interval with length, At c — 25.6msec, and in the whole 
mass parameter region, together with the parameters, t c , 
M and r\ which realize the local maximum. 



IV. TRIGGER DISTRIBUTIONS 

The data of TAMA300 contain non-stationary, non- 
Gaussian noise. Such noise cause many triggers with rate 
much larger than that expected in the stationary Gaus- 
sian noise. In order to distinguish such spurious triggers 
from triggers caused by real gravitational wave signals, 
we compute \ 2 value for each t rigg er with p > 7. The 
definition of \ 2 can be found in 17] 23] . This \ 2 is de- 
fined such that it is independent from the amplitude of 
signal if the wave form of the signal and the template are 
identical. However, since our template parameters are 
defined discretely, and thus the signals are different from 
the templates in general, when the amplitude of signal 



becomes larger, % 2 becomes larger. In order not to lose 
real signals with large x 2 , we define (, = p/ x/x 2 as a new 
statistic hj. The statistic, C, was used in our previous 
analysis [7|, and was found to be useful to distinguish 
the spurious triggers from triggers caused by real gravi- 
tational wave signals. 

The cumulative number distribution of triggers as a 
function of £ for each observation run are plotted in Fig^- 

m 



Search result: TAMA DT4, DT5 




^ threshold 

FIG. 1: The cumulative number of triggers as a function of C, 
for DT4 and DT5. 



In these plots, there are no triggers which deviate from 
the tail of the distribution of triggers significantly. This 
fact suggests that there is no candidate trigger which can 
be interpreted as real gravitational signal. 
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Search result: TAMA DT6, DT8, DT9 
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FIG. 2: The cumulative number of triggers as a function of £ 
for DT6, DT8 and DT9. 



FIG. 3: Detection probability of Galactic binaries inspirals for 
each observation run. The dashed lines show the uncertainty 
of the Monte Carlo simulations. 



V. UPPER LIMIT TO THE GALACTIC EVENT 
RATE 

In this section, we evaluate the upper limit on the rate 
of the inspirals of compact binaries in our Galaxy. In or- 
der to do this, we first evaluate the detection probability 
of Galactic signals by adding the signals to the real data, 
and by analyzing the data by the same analysis pipeline 
used in the real analysis. We assume the distribution of 
compact star binaries in our Galaxy given by 



dN = e- r2 l {2r ^e- z ' h 'rdrdZ, 



(5.1) 



where r is the radius from the center of the Galaxy, 
ro = 8.5kpc, Z is height from the galactic plane, and 
h z = lkpc. We assume that mass of each component 
star is uniformly distributed between 1 to 3 Mq , because 
we do not know much about the mass distribution model 
of binary compact stars including black holes and/or neu- 
tron stars. We also assume uniformly distributed inclina- 
tion angle of the orbital plane and the polarization angle 
of signals. The obtained detection probability is plotted 
in Fig. The DT9's data are the most sensitive to the 
Galactic events. Actually, the detection probability by 
the second half of the DT9 data is much better than that 
by the first half data. The first half data of DT9 was not 
very stable. Many triggers with large £ were produced 
by instrumental noise during that period. They degrade 
the average detection probability in DT9. 

The upper limit to the event rate from each observation 
run, Ri (i = DT6, DT8, DT9), is derived by 



Ti€i '' 



(5.2) 



where Tj is the length of data, is the detection proba- 
bility, and Ni is the upper limit to the number of event 
derived by the following formula: 



-(Ni+N- 



Elb— J! 
n=0 



e Z^n=0 nl 



1-C.L., (5.3) 



where is the observed number of triggers which ex- 
ceed a threshold, N^j is the number of triggers which are 
caused by noise alone, and C.L. is a confidence level. 

We set the false alarm rate to 1 event par year. The 
threshold which corresponds to this false alarm rate is 
evaluated by fitting the trigger distribution assuming 
that all triggers are caused by noise. We note that 
z = C 2 /2 obeys the F distribution with degree of freedom, 
(2, 2p — 2), when the data are the Gaussian noise. Here, 
p is the number of bins in the frequency region which is 
used to define \ 2 ; an d we set p = 16. In this case, the 
variable z obeys the probability density function given 
by (p — l) p (z + p— l)~ p . The cumulative number of trig- 
gers as a function of the threshold, N(z), is proportional 
to N(z) <x (p - ly-^z+p - 1)-p +1 . Thus, the plot of 
logvY(z) - ( is not linear, but logiV(z) — log(z +p — l) 
becomes linear. Although, TAMA300's data show non- 
Gaussian property, these facts suggest that log N(z) — 
log(z + p — 1) plot may be more suitable for an accu- 
rate evaluation of the false alarm rate as a function of 
the threshold. We find that this is actually the case for 
DT8 and DT9. In FigH we show the result of the fitting 
for DT9 case. The thresholds obtained in this way are 
listed in Table ITT1 On the other hand, the same plot does 
not become linear in DT6 case. We then conservatively 
select the region of the fitting so that we have a larger 
threshold for a given false alarm rate. 

With these thresholds, we have 7Y^; ; S = for all cases. 



From Eq. (|5.3() , the upper limit to the number of event is 
Ni = 2.3 for a confidence level C.L. =90%. We obtain the 
upper limit for the Galactic event rate, 130 [yr" 1 ] from 
DT6, 30 [yr- 1 ] from DT8, and 60 [yr^ 1 ] from DT9. The 
thresholds, the detection probability, and the upper limit 
for the Galactic event rate for each run are listed in Table 
ITT1 The most stringent upper limit is obtained from DT8 
data. This is because the length of data is the longest 
among three runs, and because the detection probability 
is comparable to that of DT9 on average. 

We combine these results and obtain a single upper 



6 



Search result: TAMA DT9 



E3- 



1.5 2 2.5 

Log10(z+15) threshold 



FIG. 4: The dashed line denotes the cumulative number of 
triggers as a function of log 10 (z + p — 1), (p — 16), for DT9. 
The solid line is the result of the least square fitting to this 
trigger distribution. 



limit. The upper limit from different observation runs is 
given by 



R 



iV UL 



(5.4) 



where iVuL is the upper limit to the number of events 
derived by all of data. We adopt the same threshold 
C* for each run listed in Table [Q] The total number 

of background triggers is Yli^he = 0.281. Thus, we 
have 1 event per year as a false alarm rate for combined 
DT6-DT8-DT9 data. Since the total number of triggers 
observed is zero, J2i -^obl = 0> we have TVtjl = 2.3 for 
C.L.=90%. From Eq. JO}, the combined upper limit to 
the event rate becomes 



R=17 [yr- 1 ]. 



(5.5) 





DT6 


DT8 


DT9 


Observation time [hours] 


876.6 


1100 


486.1 


Threshold (* 


21.8 


13.7 


17.7 


N (i) 


0.1000 


0.1255 


0.0555 


Detection probability 


0.18 


0.60 


0.69 


(5-Ri)fluct [yr -1 ] 


+20.6 


+2.52 


+4.04 




-24.0 


-2.82 


-3.77 


(<5-Ri)modd [yr -1 ] 


+55.4 


+4.18 


+6.84 




-16.6 


-1.53 


-2.60 




130lf 9 


JU -4.6 


60 +8.o 



TABLE II: Summary of the upper limit to the Galactic event 
rate. The errors for the upper limit are evaluated in section 
VI in detail. 



VI. STATISTICAL AND SYSTEMATIC ERRORS 

We consider various error sources which affect the de- 
tection probability. These are summarized in Table IrTTl 





1 J _L U 


DT8 


DT9 


Threshold 


+0.001 


+0.031 


+0.013 




—0.000 


-0.024 


-0.022 


Monte Carlo 


±0.093 


±0.014 


±0.080 


fall i"it" oti An 




+0.045 


+0.040 




-0.028 


-0.041 


-0.039 


(<5ei)fluct 


+0.035 


+0.056 +0.042 




-0.029 


-0.049 


-0.045 


Wave form 


-0.028 


-0.041 


-0.039 


Binary distribution 


±0.028 


±0.032 


±0.031 


model 








(<5ei)rnodcl 


+0.028 


+0.032 


+0.031 




-0.056 


-0.073 


-0.070 



TABLE III: The various error sources and their value in the 
detection probability. 

Threshold The method to derive the upper limit to the 
event rate in this paper requires the evaluation of the 
threshold which corresponds to a given false alarm rate. 
This is done by fitting the distribution of triggers as ex- 
plained in section V. There are statistical errors of the fit- 
ting due to the fluctuation of the number of background 
triggers. This error results in the error of the thresh- 
old, and the detection probability. We have the error of 
-0.02 ~ +0.03 in the detection probability in DT8 and 
DT9 cases. In DT6 explained in Section V, the 

error of the fitting due to the non-linear property of the 
distribution was already incorpolated in the fitting. The 
statistical error of the fitting itself was very small in DT6 
case, < 10~ 3 . 

Monte Carlo The error due to the Monte Carlo in- 
jeetioin test with a limited trial number is given by 
\Ai(l — £i)/ n » where rii is the number of the Monte 
Carlo trials of each run. This Monte Carlo error becomes 
about ±0.01 in the detection probability. 
Calibration The calibration of the sensitivity of 
TAMA300 is done by monitoring continuously the re- 
sponse of an injected sinusoidal test signal. The error 
of this response is much smaller than the normaliza- 
tion error discribed below, and can be neglected here. 
In the determination of the transfer function, there are 
two possible effects which affect the calibration uncer- 
tainty. One is an overall normalization error associated 
with the magnetic actuation strength uncertainty and 
its effect on calibration, and the other is uncertainty in 
the frequency-dependent response. Although the error 
in the normalization is of order 5%, the long-term drift 
is unknown. We thus conservatively adopt 10%. The 
frequency-dependent error is known to be much less than 
10%, and thus it is absorbed in the uncertainty in the nor- 
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malization. The calibration uncertainty leads to the er- 
rors of —0.03 ~ +0.05 in the detection probability. This 
calibration error is expected to depend on the different 
observation runs, and is expected to drift and/or fluctu- 
ate even within an observation runs. 
Binary distribution model We have adopted a specific 
model for the distribution of binary neutron stars in our 
Galaxy. If the distance between the Sun and the galactic 
center is different from the adopted value, the detection 
probability will be changed. The uncertainty of this dis- 
tance ±0.9kpc leads to the uncertainty of the detection 
probability about ±0.03. 

Wave form We used the wave form based on the 2.5 PN 
order. However, currently the best template has 3.5PN 
order. The uncertainty of p due to this is at most 6%. 
However, it is reported that the PN wave form itself may 
contain uncertainty [3- We thus adopt 10% reduction 
of the estimated p as an uncertainty. This produces an 
error of —0.03 ~ —0.04 in the detection probability. 

The above errors propagate to the upper limit of the 
event rate for each run. We take a quadratic sum of the 
errors due to threshold, Monte Carlo injection test, and 
calibration, since they show the property of fluctuation, 
and since they are independent each other. The sum of 
these errors, (<5ej)fl uc t, is listed in Table ITTT1 The errors 
due to the binary distribution model and the theoreti- 
cal wave form produce, simply, the shift of the detection 
probability. We thus take a linear sum of them conser- 
vatively. The sum of these errors, (8ei) mo & e \ is listed in 
Table IIIII We denote the effect of these errors to the 
upper limit for each run, Ri, as (£-Rj)a uc t and (SRi) mo< iei 
respectively, which are shown in Table ITT1 When we eval- 
uate the total error for each run, we take a quadratic 
sum of (<5i?i)fl U ct and (<5i?i) m odei, since they are indepen- 
dent each other. As shown in Table the errors of 
the upper limit to the event rate for each run become 
+59/ - 29 [yr- 1 ] for DT6, +4.9/ - 4.6 [yr" 1 ] for DT8, 
and +8.0/ - 4.6 [yr" 1 ] for DT9. 

Finally, we evaluate the error for the combined upper 
limit, Ea. l|5.4|l . The effect of (Sei) to R is evaluated by 
taking a quadratic sum of each effect of (Sei) to R, and 
we have +0.965/ — 1.08 [yr -1 ]. The effect of (<5ei) m odci to 
R is evaluated by simply shifting each ej in Eq. I|5.4|) , and 
we obtain +2.86/ — 1.05 [yr -1 ]. The total error in R is 
evaluated by taking a quadratic sum of these two errors. 
We have the upper limit with error, R = 17l'J'g^[yr -1 ]. 
By taking larger value as a conservative upper limit, we 
obtain 

i? = 20[yr- 1 ]. (6.1) 

This value is much larger than an astrophysically ex- 
pected value, 8.3 x 10 -5 [yr _1 ] for the coalescence of 
neutron star binaries. However, this rate is smaller than 
that obtained by LIGO S2 search, 47 [yr^MWEG -1 ], 
or by LIGO-TAMA joint analysis, 49 [yr^MWEG -1 ]. 
Main reason for this is that the length of data used in 
our analysis is much longer than these analyses. 



VII. SUMMARY AND DISCUSSION 



In this paper, we have presented the results from 
the TAMA300 data analysis to search for gravitational 
waves from inspiraling compact binaries in a mass range, 
1 - 3M . We analyzed DT4, DT5, DT6, DT8, and DT9 
data of TAMA300. There were no triggers which deviate 
from the tail of the distribution of triggers significantly. 
We thus conclude that there is no candidate trigger which 
can be interpreted as a real gravitational signal. By us- 
ing the long and sensitive data from DT6, DT8 and DT9, 
we obtained upper limits to the Galactic event rate from 
each observation run. We combined these results and 
obtained a single upper limit, 20 [yr" 1 ] at a 90% confi- 
dence level from these three observation runs. We evalu- 
ated the systematic errors due to various effects such like 
the uncertainty of calibration and the uncertainty of the 
background estimation. In the upper limit, these effects 
are included. 

The upper limit obtained in this paper is much larger 
than an astrophysically expected value for the coales- 
cence of neutron star binaries. However, this upper limit 
is significant since it is derived by observation. Never- 
theless, more sensitive detectors are necessary to obtain 
more stringent upper limit to the event rate, and to de- 
tect the signal. TAMA300 is now improving the suspen- 
sion system by installing the Seismic Attenuation System 
in order to obtain better sensitivity and better stability. 
When it is finished, it is expected to have much better 
sensitivity than DT9. LIGO has already been performed 
3rd and 4th scientific runs with better sensitivity than 
S2. Further, LIGO is now conducting the 5th scientific 
run since November 2005, with its design sensitivity. It 
can detect the inspiraling binaries up to ^lOMpc dis- 
tance. They are expected to be able to set a much more 
stringent upper limit. 

When the spin angular momentum of compact objects 
cannot be neglected, the spinless template is not good 
enough to detect the signal, and we need to employ tem- 
plates with spins. However, since the number of parame- 
ters becomes much larger than the non-spinning case, it 
requires very powerful computer resources. One way to 
avoid the use of the full templates with spins will be to 
use some phenomenological templates with small num- 
ber of parameters [2(j ■ We will work on such cases in the 
future. 

Despite of the improvement and long term observation 
of current detectors, the chance to detect gravitational 
waves by these first generation detectors will not be very 
large. We need more sensitive detectors, such like ad- 
vanced LIGO UK and LCGT ^ . These detectors will 
detect gravitational waves frequently, and will be used 
to investigate the strong field region of gravity and the 
astrophysics of compact objects. 
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